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Particle velocity profiles upon shock compression and isentro- 
pic release were measured for polycrystalline calcite (Solenhofen 
limestone) to 12 to 24 GPa and for porous calcite (Dover chalk, 
p 0 = 1.40 g/ cm 3 , 49% porosity) to between 5 and 11 GPa. The elec- 
tromagnetic particle velocity gauge method was used. Upon shock 
compression of Solenhofen limestone, the Hugoniot y, as tic limit 
was determined to vary from .36 to .45 GPa. Transition shocks at 
between 2.5 and 3.7 GPa, possibly arising from the calcite II-III 
transition, were observed. For the Solenhofen limestone, the 
release paths lie relatively cIosk to the Hugoniot, Evidence for the 
occurrence of the calcite III-I1 transition upon release was 
observed, but no rarefaction shocks were detected. Initial release 
wave speeds suggest retention of shear strength up to at least 20 
GPa, with a possible loss of shear strength at higher pressures. The 
measured equation of state is used to predict the fraction of 
material devolatilized upon isentropic release as a function of 
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shock pressure. The effect of ambient partial pressure of C0 2 on 
the calculations is demonstrated and should be taken into account 
in models of atmospheric evolution by means of impact-induced 
mineral devolatilization. Mass fractions of C0 2 released expected 
on the basis of a continuum model are much lower than deter- 
mined experimentally. This discrepancy, and radiative characteris- 
tics of shocked caicite indicate that localization of thermal energy 
(shear banding) occurs under shock compression. Release isen- 
trope data indicates that Dover chalk loses its shear strength when 
shocked to 10 GPa pressure. At 5 GPa the present data are ambi- 
guous regarding possible shear strength. For Dover chalk, shock 
entropy calculations result in a minimum estimate of 90% devola- 
tilization upon complete release from 10 GPa. Isentropic release 
paths from calculated continuum Hugoniot temperatures cross 
into the CaO (solid) + C0 2 (vapor) field at improbably low pressures 
(for example 10~ 7 GPa for a shock pressure of 25 GPa). However, 
release paths from measured shock temperatures cross into the 
melt plus vapor field at pressures greater than .5 GPa, suggesting 
that devolatilization is initiated at the shear banding sites. 
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Introduction 

The shock compression and release behavior of carbonates is of interest in 
the study of cratering mechanics, shock metamorphism in carbonate terranes 
and the generation of C0 2 bearing atmospheres on the terrestrial planets. 
Approximately 30% of the known or probable terrestrial meteorite impact 
craters occur at least partially in carbonate rocks [Grieve and Robertson, 1979]. 
Impact-induced devolatilization of hydrous and carbonate minerals appears to 
play a role in the evolution of terrestrial planetary atmospheres [Lange and 
Ahrens, 1983, 19t'4]. 

Static compression studies indicate that two metastable polymorphs of cal- 
cite (CaC0 3 ), exist at pressures above that of the calcite-aragonite transition 
(Bridgman. 1939; Jamieson, 1957). The single crystal calcite Hugoniot shows 
density discontinuities at pressures of between 1.8 to 2.4 GPa, at about 3 GPa, at 
about 4.5 GPa, and at about 9.5 GPa, corresponding to the Hugoniot elastic limit 
(HEL) and the calcite I-Il, II-1II, and 1II-V1 (see below) transitions, respectively 
[Ahrens and Gregson, 1964]. Analogous transitions occur in shocked limestone 
(polycrystaiiine calcite), although at somewhat lower pressures [Ahrens and 
Gregson, 1964, Grady et a l., 1978]. The aragonite Hugoniot also displays evi- 
dence of several transitions [Vizgirda and Ahrens, 1982], which are probably 
unrelated to calcite II and III since aragonite is denser than these phases. Above 
about 10 GPa the single crystal calcite [Ahrens and Gregson, 1964], polycrystal- 
line calcite [Adadurov et ai., 1961; Kalashnikov et a L, 1973; van Thiel et a t., 
1977], and single crystal aragonite [Vizgirda and Ahrens, 1982] Hugoniots are 
very similar, suggesting transformation to a similar high pressure polymorph 
above that pressure. In the remainder of the text we refer to this stable high 
pressure polymorph as calcite VI. Calcite IV and V are low-pessure, high- 
temperature forms of CaC0 3 [Carlson, 1983]. 
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Release isentropes for polycrystalline calciie rocks shocked to pressures up 
to about 4 GPa have been determined using laser interferometry [Schuler and 
Grady, 1977; Grady at ed., 197B; Grady. 1983] and electromagnetic particle velo- 
city gauges [Larson and Anderson. 1979]. Murri at a l. [1975] employed elec- 
tromagnetic particle velocity gauges to determine release paths for selected 
carbonate rocks (porosity 0 and 15%) shocked to between 10 and 30 GPa. In this 
pressure range the initial release paths lie at greater densities than the 
Hugoniot. Release paths for single crystal aragonite shocked to pressures up to 
40 GPa have been determined using buffer and inclined mirror techniques [Viz- 
girda and Ahrens, 1982]. For shock pressures up to about 13 GPa the reirase 
pat__s are steep, ana maximum post-shock densities are greater than the initial 
volume. For shock pressures, between about 13 and 40 GPa the release paths 
generally lie close to the Hugoniot, and maximum post-shock densities are less 
than or equal to the initial density. 

Shock loaded calcite decarbonates (i.e., decomposes to CaO plus C0 2 gas) 
when shocked to pressures as low as 10 GPa [Lange and Ahrens, 1983, 1984]. 
However, calculations of post-shock energy content assuming that the release 
path is identical to the Hugoniot require a shock pressure of 45 GPa for incipient 
devolatilization [Kieffer and Simords, I960], and consideration of post-shock 
entropy content assuming isentropic release [Zel’dovich and Raizer, 1967; 
Ahrens and O’Keefe, 1972] requires shock pressures of 33 GPa for incipient devo- 
latilization [Vizgirda and Ahrens, 1982]. However, these previous applications of 
post-shock energy and entropy calculations have failed to consider the effect of 
the ambient partial pressure of the volatile species on the equilibrium of the 
devolatilization reaction. 

A shock temperature measurement at 40 GPa on single crystal calcite 
yields a color temperature of 3700 K, which is over twice the Hugoniot 
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temperature calculated from continuum thermodynamic models, with an emis- 
sivity of 0.0025 [Hondo and Ahrens, 1983]. This result has been interpreted as 
indicating the presence of a large number of closely spaced high temperature 
shear-band regions immediately behind the shock front, in support of the shear 
instability models of Grady [1977, 1980] and Horie [1980]. These zones of 
intense heating would cause devolatilization to begin at much lower pressures 
than those calculated on the basis of continuum models. 

In an effort to determine the limits of applicability of continuum models, we 
have determined the release paths for a slightly porous polycrystalline car- 
bonate (Soienhofen limestone, p 0 = 2.58 g/ cm 3 ) shocked to pressures between 
12 and 24 GPa and for a very porous calcite (Dover chalk, p 0 = 1.40 g/ cm 3 ) 
shocked to between 5 and 11 GPa pressure using the electromagnetic particle 
velocity gauge technique. Release paths have not been previously determined 
for these materials shocked to these pressures. The release path data are used 
to determine the post-shock energy gain which is then compared to estimated 
values of energy input required for incipient volatilization Ejy. In addition, the 
sound velocity data obtained for the shocked state and along the unloading path 
provide constraints on the mechanical strength properties of these materials 
during shock compression and unloading. 

Experimental Technique 

The particle velocity gauge technique used in these experiments is 
described elsewhere in detail [Boslough, 1983], and shown schematically in Fig- 
ure 1. The targets consisted of four 1.5 mm thick plates of Soienhofen limestone 
with 12.5 fim thick copper, polyamide-backed (Kapton, 12.5 fxm) particle velo- 
city gauges epoxied between each plate and onto the free surface. The faces of 
the assembled target were parallel to within 25 fx m. Archimedean and bulk den- 
sities were determined individually for each plate. The chalk target assemblies 


- 8 - 


f 


are shown in Figure 2. In the chalk assemblies the gauges were stretched across 
each plate and epoxied only on a 3 mm wide strip along the plate edge. Because 
the chalk is relatively weak, each target assembly consisted of only three 2 to 
2.5 mm thick plates. Additional clamping was provided by a circular clamp ring 
on the front and clamping strips on the rear of the mounting ring. 

The target was mounted at the center of a set of Helmholtz coils so that, the 
active element of each gauge, the incident projectile velocity, and the magnetic 
field generated by the coils were mutually perpendicular. Upon impact a voltage 
is induced across the active element of the gauge which is proportional to the 
effective gauge (active element) length L. the magnetic induction B. and the par- 
ticle velocity u p . 

V(t) = BLu p (t) (1) 

where t is the time after impact. V(t) is in units of volts when B is in gauss, L is in 
cm. and Up is in cm/sec. In these experiments B was approximately 1.8 kG and 
L was approximately 0.95 cm. The gauge voltages were recorded by a series of 4. 
Tektronix oscilloscopes using the circuits sketched in Fig. 3. 

The shots were performed on the 40 mm propellant gun at the California 
Institute of Technology [Ahrens ef aZ.. 1971]. The impact tank pressure was 
between 50 and 160 fxm Hg. Fused quartz and polycrystalline alumina flyer 
plates between 6 and 18 mm thick were employed. The projectile velocity was 
measured by determining the flight distance of the projectile during the known 
time interval between two flash x-ray photographs of the projectile. 

Typical experimental results are shown in Figure 4, which shows particle 
velocity versus time after impact obtained from digitized oscilloscope records. 
The shock wave generated by the impact reaches each of the gauges in sequence 
(gauge 1 is nearest the impacted surface, gauge 4 is on the rear, or free, 
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surface) accelerating it to the particle velocity of the Hugoniot state. Gauge 4 is 
accelerated immediately to the free surface velocity. Gauge 3 and then gauge 2 
are later accelerated to higher particle velocities by the rarefaction wave pro- 
pagating back into the sample from the free surface. Gauge 2 is disrupted by 
rarefactions originating at the edges of the target before reaching the final free 
surface velocity. Gauge 1 is intercepted by the forward-travelling rarefaction 
originating at the upstream side of the flyer plate and is therefore accelerated 
to lower particle velocities. Theie relationships are shown more clearly in 
schematic position-time plot, shown in Figure 5. 

The time interval between the shock arrival at each gauge was used to 
determine the shock wave velocity. This was combined with the projectile velo- 
city and the known flyer material Hugoniot [Marsh, ll>80] in an impedance 
match solution [McQueen at a L. 1970] to yield the Hugoniot state. The Hugoniot 
particle velocity obtained in this way provided the reference point for the gauge 
voltage versus particle velocity calibration. 

The stress and density along the release path are given by equations for 
conservation of mass and linear momentum [Cowperthwaite and Williams, 1971] 

V' a "A = £-- c “ <2> 

( da/ d u? ] h = p 0 C„ (3) 

where p is the density, p 0 is initial density, Up is particle velocity, h is the 
Lagrangian space coordinate, and 

C u = (ah/ dt) u (4a) 


c„ = (ah/ at). 


(4b) 
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where t is time. C„ is the velocity of propagation of a wave with particle velocity 
Up. and C ff is the velocity of propagation of a wave with stress o. The release 
waves are nonsteady simple waves [Courant and Friedrichs, 1948; Petersen at 
at.. 1970} Therefore C u and C e are equal, but depend on Up. The digitized parti- 
cle velocity as a function of time records were used to compute C u and then 
equations (2) and (3) were integrated numerically to obtain the stress-density 
path of the release. The Lagrangian sound speed, C u . is obtained using the finite 
difference approximation, 


C„ 


Ah 

At 


( 5 ) 


where Ah is the initial distance between gauges and At is the transit time for a 
disturbance with a particle velocity u p . Eulerian sound speeds are equal to 

C E = ^-C u (8) 

The Eulerian sound speed corresponds to the sound speed relative to the labora- 
tory reference frame. 

Results 

The experimental results for the shocked states of the Solenhofen lime- 
stone and the Dover chalk are summarized in Table 1. Figures 6 and 7 are shock 
velocity-particle velocity (U,-Up) and pressure-density (P-p) plots of the data, 
respectively. The experimental data for shot 599 (Dover chalk) appear to be 
erroneous. The results are listed in Table 1, but are not plotted. No release 
path data were obtained from shot 601. 

In two of the shots, shots 598 and 602 (Hugoniot pressures 13.10 and 11.96 
GPa, respectively), multiple wave structure was observed (Figure 8). Transitions 
occurring between .36 and .45 GPa and between 2.5 and 3.7 GPa were identified. 
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Multiple wave structure has be .n observed in previous shock wave experiments 
on carbonates, and is caused by dynamic yielding and at least 3 additional phase 
changes that calcite undergoes at low pressures [Ahrens and Gregson, 1984; 
Grady ef a £., 1978]. However, in the experiments reported here the shock and 
particle velocities of the intermediate states are poorly constrained, and are not 
completely consistent. These less reliable data are shown in parentheses in 
Table 1. 

The velocity of 5.7 ± 0.3 km/sec for the first wave of shot 802 (Table 1) is 
comparable to the longitudinal wave velocity of Solenhofen limestone, 5.84 - 5.97 
km/sec [Peselnick, 1982; Hughes and Cross, 1951] indicating that this transition 
represents the Hugoniot elastic limit (HEL). Ahrens and Gregson [1964] report 
elastic precursor velocities of 5.3 to 5.7 km/sec and a HEL of 1.0 to 1.5 GPa for 
carbonate-bearing rocks. Grady ef a l. [1970] report a break in the Solenhofen 
limestone loading wave profile at 0.8 GPa and conclude that the calcite Ml phase 
transition occurs at this pressure, coincident with the onset of dynamic yielding. 
Our data are not of sufficient resolution in this pressure range to resolve this 
question. Comparison with the results of Grady et a L. [1978] suggests that the 
transition at 2.5 to 3.7 GPa corresponds to the calcite II-1II transition. In the 
higher shock pressure shots these transitions are apparently overdriven, in 
agreement with previous shock experimental results [Ahrens and Gregson, 
1964]. No precursory waves were observed in the shots on Dover chalk. 

Figure 6 shows the U,-Up data determined in this study compared to trends 
observed in various calcium carbonate -bearing rocks and minerals by previous 
investigators. For Solenhofen limestone the values of U, and u p lie within the 
range of values previously determined by other investigators, but the slope of 
the line defined by the data is greater. The data in the range of 12 to 24 GPa can 
be fit by an equation of the forir 1 
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u, = c 0 + a Up (7) 

where C 0 and s are constants, with C 0 = 3.289 km/sec and s = 1.798 (r* = 0.97). 
Previously determined values of C e and s range from 3.62 to 3.99 km/sec and 
1.32 to 1.61, respectively, and are summarized in Table 2. The U,-Up data for 
Dover chalk, p 0 = 1.40 g/cm 3 , can be fit with C 0 = 0.687 and s = 1.598 (r 2 = 
0.989). This line is consistent with the trend formed by the porous calcite 
results of Kalashnikov et a l,, [1973]. 

As shown in Figures 4 and 8, the particle velocity records of all the Solenho- 
fen limestone shots indicate that a phase transition occurs upon release. The 
pressures and densities at which it occurs for each shot are listed in Table 1. 
The data are somewhat scattered but the transition occurs at about 5.5±2.5 
GPa, and no trend with peak stress level is apparent. Previous observations of a 
phase transition upon release in this pressure range have been attributed to the 
calcite Ill-II transition [Grady et a L, 1978]. Rarefaction shocks have been 
observed in carbonates shocked to these pressures [Murri et a l., 1975]. How- 
ever. no evidence for rarefaction shocks was seen in these experimental 
records. 

The initial release paths for Solenhofen limestone plotted in Figure 7 lie 
near the Hugoniot, but at slightly higher densities, for all shock pressures stu- 
died. The particle velocity records do not give complete information down to 
zero pressure due to disruption of the foil gauges by rarefactions propagating 
inward from the sample edges. Extrapolation to zero-pressure yields complete 
release densities between 2.59 and 2.81 g/cm 9 . In contrast to the behavior for 
Solenhofen limestone, release paths for single crystal aragonite vary 
significantly with shock pressure. For shock pressures between about 10 and 15 
GPa the aragonite release paths are very steep, with complete release densities 
(determined via inclined mirror experiments) equal to or greater than the initial 
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density. For shock pressure greater than about 20 GPa, the release paths 
approximate the Hugoniot, with complete release densities equal to or less than 
the initial density [Vizgirda and Ahrens, 1982]. Release paths for Solenhofen 
limestone [Schuler and Grady, 1977] and other carbonate rocks [Grady at a 
1978; Grady, 1983] shocked to pressures up to about 5 GPa generally follow the 
Hugoniot, with extrapolated complete release densities approximately equal to 
the initial density. 

The release paths for Dover chalk lie at higher dersities than the Hugoniot. 
The extrapolated complete release densities lie between 2.44 and 2.55 g/c m 3 . 

Eulerian sound speeds for the hocked states were calculated according to 
equations 5 and 8 and are plotted in Figure 9. 

Discussion 

Properties of calcite VI. reached at about 10 GPa shock pressure, have been 
deduced by Vizgirda and Ahrens [1982] and are listed in Table 3. The zero- 
pressure bulk sound velocity 


V, = (Ko^p 0 ) 1/8 


(8) 


where K*, is the zero pressure isentropic bulk modulus, is between about 5.00 
and 5.54 km/sec at 300 K. Assuming a Poisson’s ratio a of 0.25 we obtain a 
zero-pressure, 300 K compressional velocity V p of between 8.48 and 7.15 km/sec 
using the relationship 


V p s 


K 0# /p-[l J -2 (l-2o)/ (1 + p)] 


i/z 


(9) 


We can compare the observed Hugoniot sound velocities with compressional 
wave velocities V p along an isentrope for calcite VI using finite strain-theory 
[Sammis at a l.. 1970; see also Burdick and Anderson. 1975]. Thus. 
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vj (p) * v*. (1-2*) [l-2*(3 D, - 1)] (10) 

V* (p) * V«, (l-2e)[l-2e(3 D. - 1)1 (11) 

and 

P * -3K* e (1-2*) 8 ' 8 (1+2* ()• (12) 

The volumetric strain e is 

e = [l-(p/p 0 ) 8/8 ]/2 (13) 

the finite strain parameter ( is 

( = 3[4-(dK 0 ^dP),]/4 (14) 

D, = Ko.dinVp/dP (15) 

and 

D.= K^dinV./dP (18) 

No determinations of D p and D, for calcite VI exist. For many minerJs D p lies in 
the range 1.0 at D p « 1.6 and D ( lies between 0.1 £ D, < 1.0 [Anderson ef a l , 1968; 

Sammis *t a l., 1970]. This range of values was used in computing 

Vp(p) and V, Ip). The variation of bulk sound speed V, along an isentrope is given 
by 

V 8 (p) =Vp(p)«-4V*(p)/3 (17) 

The effect of temperature on sound velocity is small. The temperature 
along the Hugoniot T H was calculated using continuum methods [Zel'dovich and 
Raizer, 1967] and is listed for each shot in Table 4. The parameters used for the 
calculation of Th at e listed in Table 3. For the Dover chalk, the temperatures 
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are sufficiently high that calculated Hugoniot temperatures are relatively insen- 
sitive to the specific isentrope employed in the calculation. Assuming a value of 
(d Vp/ dT) p of -3.3 x 10~ 4 1cm/ sec * K (values tabulated by Anderson ef of. [i960] 
range from about -1.3 x 10“ 4 to -3.2 x 10~ 4 km/ sec • K), an increase in the initial 
temperature of 430 K results in a Vp decrease of 0. 13 km/sec. 

The range of values of V p and V, along a 730 K isentrope calculated using the 
parameters listed in Table 3 is shown in Figure 9. Although the calculated bulk 
and compressions! sound speeds are poorly constrained, the data indicate that 
the measured initial release wave speed is equal to that of the compressions! 
wave velocity for shock pressures up to about 20 GPa. In other words, calcite VI 
retains its shear strength in the shocked state up to this pressure. At the 
highest pressures attained. 22 to 24 GPa, a loss of strength in the shocked state 
may occur. 

The bulk sound speed V # at the Hugoniot may also be emulated from the 
shock wave equation of state (equation (7)) and a model for the Gruneisen 
parameter y, using the relationship 


v,= 


= V H I — I 
H dV 1 


hi 


(Vc-V H ) 



1/2 


♦Pk7/2V H 


( 10 ) 


[.McQueen *t al „ 1987] where V e is the initial sample volume, V H is the Hugoniot 
dP 

volume, and i* the slope of the Hugoniot. We use the relationship 

7 s 7o (p©/p) w here y 0 * 1.5 and n * 1.0 [Virgirda k Ahrens, 1982] (Table 3). This 
curve is plotted as a dashed line in Figure 9 and supports the conclusion that a 
loss of shear strength occurs at a pressure of 22 to 24 GPa. 

The Eulerian sound speeds along the release paths of the shocked Dover 
chalk are plotted in Figure 10. At the relatively low shock pressures and high 
continuum temperatures reached in these experiments it is unlikely that the 
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starting material transformed to calcite VI. However, it is not known which of 
the three calcite low pressure polymorphs exists in the shocked state. We have 
therefore plotted the bulk and compressional wave velocities for calcite 1, II, and 
III using, for comparison, the data of Singh and Kennedy [1974], assuming that 
Poisson’s ratio a - 0.25 [Anderson si a l., 1988; Grady ef oZ., 1978]. Also plotted 
is the bulk sound speed in the shocked state calculated using the low pressure 
(5-13 GPa) polycrystalline calcite Hugoniot of Adadurov et al. [1961] and equa- 
tion (18). The datum for shot 620 plots between the bulk and compressional 
wave velocities for calcite I and III, whereas the result for shot 622 (Ph = 9.47 
GPa) lies near the bulk sound speed of calcite I and III. 

The behavior observed for the limestone and chalk may be compared with 
observations on other minerals. San Gabriel anorthosite shocked to pressures 
between 6 and 10 GPa retains its shear strength in the shocked state, but under- 
goes a strength loss upon release [Boslough and Ahrens, 1984]. Quartz shocked 
into the mixed-phase region [Grady at al., 1975] and non-porous carbonates at 
shock pressures up to about 5 GPa [Grady et al., 1978] have initial release wave 
speeds equal to the bulk sound speed. There appear to be two mechanisms that 
cause loss of shear strength in the shocked rtate. The onset of shear strength 
loss occurs upon the initiation of phase transition. Thus, materials undergoing 
phase transitions along the Hugoniot (quartz, carbonates at low pressure) exhi- 
bit initial release wave velocities equal to bulk sound speed, whereas materials 
undergoing transformation upon release (carbonates above 10 GPa shock pres- 
sure, anorthosite) exhibit initial release wave velocities corresponding to 
compressional wave speed, with a reduction to lower velocities upon release. 
Shear strength loss also accompanies extensive shear band formation [Grady et 
al., 1975], A‘ 22 GPa Solenhofen limestone exhibits a loss of shear strength. At 
this shock pressure, the measured shock temperature is approximately 2455 K 
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(D. Schmitt, personal communication), over twice the calculated continuum 
temperature. Apparently, at this shock pressure the shear band density is large 
enough to affect the bulk material properties. 

Shock recovery experiments on single crystal calcite demonstrate that the 
onset of decarbonation according to the reaction 

CaCOjnr. impact Ca0 («oM) + C°2 (im) (19) 

occurs upon release from shock pressures of less than 10 GPa (Lange and 
Ahrens, 1983). Jpon complete release from 12-13 GPa shock pressure, about 30 
per cent of the C0 2 initially present is volatilized, and 50-60 per cent of the C0 2 
is lost upon release from 20-25 GPa [Lan^. Sc Ahrens, 1983]. In the following sec- 
tions we compare the experimentally determined values of the fraction of C0 2 
volatilised with calculations based on post-shock energy and entropy content. 

The post-shock energy increase AE p in a material shocked to Hugoniot pres- 
sure Ph and volume Vr is given by [Ahrens and O’Keefe. 1972, Zel’dovich and 
Raizer, 1967] 

V oo 

ASp = VEo = |"(P. + Ph) (V. - V H ) - / (PdV). (20) 

’H 

where Ep is the post-shock energy content, E e is the initial energy content, V 0 is 
the initial specific volume, and V 00 is the complete release volume. The integral 
is performed along the release path, which is assumed to be isentropic [Kieffer 
and Delaney, 1979; Jeanlcz and Ahrens, 1979]. The first term on the right hand 
side of equation (20) is the internal energy increase upon loading the sample to 
the Hugoniot. and the second term represents the work performed in expanding 
the material to V 00 along the measured Please path. The experimental release 
path results do not extend to zero pressure, so V 00 was estimated by computing 
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an isentrope from the end of the measured release path to zero pressure. For 
the Limestone shots the calcite VI parameters were used; for the chalk shots cal- 
cite I and II parameters were used (Table 3). This approximation leads to an 
underestimate of the integral term in equation (20), and therefore to an overes- 
timate of &E P . The assumed complete release densities and calculated post- 
shock energy increases are listed in Table 4. Figure 11 is a plot of post-shock 
energy increase AE p as a function of shock pressure. Although the data are 
somewhat scattered, a trend toward higher post-shock energy content with 
increased pressure is observed. For the Solenhofen limestone AE P reaches 300- 
400 J/g at about 25 GPa. For the Dover chalk shocked to 10 GPa, AEp is about 
1600 J/g. The curves in Figure 11 were calculated assuming isentropic release 
from the Hugoniot. For the Solenhofen limestone shots, the zero-pressure, isen- 
tropic release temperature T,*] was calculated using the parameters for calcite 
VI, whereas for the Dover chalk shots the calculation is based on the calcite 1 
and calcite II parameters (Table 3). Because of the large amount of shock heat- 
ing occurring upon compression of the initially porous chalk, post-shock energy 
auid entropy values are only moderately sensitive to the choice of isentrope. The 
calculated post-shock energy gain AE p was then obtained using 

^r*l 

AEp = £ C v d T (21) 

where T 0 is the initial temperature and C v is the calcite heat capacity at con- 
stant volume. The integration is carried out along the 1 bar isobar. C v was calcu- 
lated using the C p (heat capacity at constant pressure) values of Robie et al. 
[1970]. The calculated values of post-shock energy content are in good agree- 
ment with the experimentally determined values. 

The energy required for incipient decarbonation (volatilization) Ejy can be 
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calculated using the data of Robie at al. [1978]. It is important to bear in mind 
that the equilibrium of the reaction 

CaC0 3 (eiicit*) - CaOjioiy) + C0 2 [gtm) (22) 

is strongly influenced by the ambient pressure of C0 Z , P C0( . Thus for pure solid 
calcite and pure solid CaO, the Gibbs free energy of reaction (22) A G n is 

A = A G£ + RT£n (23) 

where the reference state for the solids is defined as pure solids at 1 bar pres- 
sure and the temperature of interest T, and for C0 2 the reference state is pure 
C0 2 gas at 1 bar (10 s Pa) pressure and T. Thus for Pco 8 - 1 bar the temperature 
of incipient vaporization Tjy is approximately 1171 K and the energy of incipient 
vaporization Ejy, given by 

T iv 

EjV = ( ca jcit«) dT (24) 

is approximately 980 J/g [Vizgirda and Ahrens, 1982; Kiefler and Simonds, 1980]. 
However, for P^ « 3.3xl0 -4 bar (corresponding to dry air at 1 bar total pres- 
sure ) Tiv is approximately 800 K, and Enr is about 530 J/g, and for Pco g ~4.3xl0~ fl 
bar (corresponding to dry air at a total presure of 100 fim Hg) T^ is about 596 K 
and Ejv is 296 J/g. Thus Eiv is strongly dependent on the explicit value of Pco 8 
used in the calculation of Tiv, and conversely, the amount of C0 2 evolved during 
release may depend on the ambient P COe . 

Such equilibrium considerations do not apply absolutely in a dynamic shock 
experiment, but evidence that ambient C0 2 pressure affects the amount of 
evolved C0 2 comes from comparison of the work of Boslough et al. [1982] with 
that of Lange and Ahrens [1983] and that of Kotra et al. [1983]. The former 
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investigators studied the evolution of C0 2 during shock compression of calcite by 
capturing the evolved gas in an initially evacuated, confined chamber. For a 
shock pressure of 16 GPa only 0.03 to 0.3 per cent of the C0 2 was volatilized, and 
the final d C q 1 in the chamber at room temperature was between about 2.7xl0~ 3 

to 7.2xl0 -3 bars [M. Boslough, personal communication, 1984]. For Pco 2 = BxlO -3 
bars. T iv = 900 K and £{v - 648 J/ g. The studies of Lange and Ahrens [1963, 
1984] were performed using vented assemblies [Lange and Ahrens. 1982]; there- 
fore the final Pco g was that of air. 3.4xl0 -4 bar. For 
p co = 3.4X10" 4 bar. Tjy = 800 K and Ejv = 528 J/g. They found significant decar- 

w 

bonation; as much as 30 per cent at pressures as low as 10 GPa (see Figure 12). 
The analytical method was different in each of the two sets of experiments; 
Boslough et a l. determined the amount of C0 2 evolved whereas Lange and Ahrens 
determined the amount of C0 2 remaining in the recovered solid material. 
Nonetheless, it appears that at least part of the difference in results may be 
attributed to differences in P COg . The results of Kotra ef al [1983] differ from 
those discussed above because their targets (both vented and non-vented) were 
exposed to the COg-rich muzzle gases from the gun barrel. The very low fraction 
of C0g volatilized as a function of shock pressure obtained by these workers may 
be caused by very high Pco g during release. 

There are several implications of these results. 

1. The release paths determined in this study may not be completely com- 
parable to the volatilization results of either Boslough et al. or Lange and 
Ahrens, inasmuch as the experiments reported here were performed at about 
100 /*m air pressure (Pco g ~4.3xl0*°bar). 

2. The evolution of Pcog must be taken into account when calculating the 
evolution of impact-induced planetary atmospheric C0 2 content. Note also that 
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these considerations apply to the shock entropy criterion for vaporization upon 
release [Ahrens and O'Keefe, 1972; Zel’dovich and Raizer, 1967], because this 
method also depends on the identification of a Tiv for calculation of Sjy (entropy 
of i :ipient vaporization). Furthermore, the equilibrium considerations apply to 
any shock devolatilization process, especially shock dehydration [c./. Boslough 
et a £., 1980; Lange and Ahrens, 1982]. 

The s*-3ek entropy criterion for incipient or complete vaporization 
[Zel’dovich and Raizer, 1967; Ahrens and O'Keefe, 1972] is convenient to use 
because it does not require experimental knowledge of the exact release path. 
Instead, the entropy gain in the shocked state ASh is calculated according to 

AS H = S tr +C,in(T H /T,) (25) 

where is the entropy of transition from calcite to calcite VI at 1 bar and 298 K 
and T h and T, are calculated by continuum methods (Zel'dovich and Raizer, 
1967). The entropy of transition has been estimated to be between -.116 and 
-.086 J/g on the basis of entropy-molar volume systematics [Vizgirda and 
Ahrens, 1982] and we have adopted the value of -0.101 J/g for use here (Table 3). 
The shock entropy is dependent on T H , which in turn is greatly influenced by the 
initial porosity of the material. Assuming isentropic release from the shocked 
state [Kieffer and Delaney, 1979; Jeanioz and Ahrens. 1979] the entropy of the 
completely released state will be equal to that of the shocked state, and this 
value is compared to the entropy of incipient vaporization Sjy. where 

T iv C 

Srv = ^ -^-dT (26) 

where C p is the atmospheric pressure heat capacity at constant pressure, 
obtained from Robie ef ai. [1978]. As discussed above, Tjy will be dependent on 



- 20 - 


th e P C0| chosen. The entropy of complete vaporization at Tjy, S ev . is computed 
from the values of entropy as a function of temperature in Robie at al. [1976]. 
Figure 13 is a plot of shock entropy versus shock pressure for calcite of varying 
porosities calculated using the calcite Hugoniot of Adadurov at al. [1961] and the 
parameters for calcite VI listed in Table 3. The Hugoniots for the porous 
material are calculated according to 


1-2. liS--! 
2 V H 


l-Z-IftL-i 

2 V H 


where P p and Pd are the pressures at volume V« of the porous material and the 
dense material, respectively, y is the Gruneisen parameter, and V CiP and V 0 jq are 
the zero-pressure volumes of the porous and dense materials, respectively 
(Zel’dovich and Raizer, 1967). The experimentally determined shock entropies 
from this study are also plotted. The results for the Solenhofen limestone (ini- 
tial porosity about .04) indicate that Sjy for Pco n = 4.3xl0 -fl bar (dry air atmo- 

• 

sphere at 100 fjm Hg total pressure, the conditions of the experiment) is 
reached at about 12-20 GPa, and that at about 25 GPa between 20-40% of the C0 2 
would be volatilized. The results from the Dover chalk shots (initial porosity of 
0.49) indicate that greater than 90% of the material would be decarbonated at 10 
GPa shock pressure under the experimental conditions. For an air atmosphere 
(Pcoj = 3.3xl0 -4 bar) the figure indicates a minimum pressure for incipient devo- 
latilization of non-porous calcite of 21 to 32 GPa, compared with the experimen- 
tal value of less than 10 GPa determined by Lange and Ahrens. 

The shock energy or entropy criteria for incipient and complete vaporiza- 
tion when applied using the ambient (or final) vapor composition yield the lowest 
values of Ejv or Sjy and thus predict the greatest amount of devolatilization con- 
sistent with equilibrium thermodynamics. Thus, equilibrium thermodynamic 
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considerations fail to describe the experimental results for Solenhofen lime- 
stone even when the most generous assumptions concerning the final state of 
equilibrium are made. 

We therefore conclude that a mechanism for the localization of thermal 
energy such as the shear instability model [Grady. 1977, 1980; Horie, 1980] is 
required to explain the experimental observations. In the case of single crystal 
calcite the existing shock temperature determinations support this conclusion. 
Kondo and Ahrens [1983] report a shock temperature of 3700K and an emissivity 
of 0.0025 for single crystal calcite shocked to 40 GPa. This temperature is 
several times the continuum temperature of 1300-1500K [Vizgirda and Ahrens, 
1982]. Textural examination of recovered calcite shocked to 40 GPa suggests 
the presence of partially molten material into which shock released C0 2 has 
been injected [Lange and Ahrens, 1984]. Schmitt (personal communication) has 
measured a shock temperature of 2455 K for single crystal calcite shocked to 
22.5 GPa. 

The relationships between measured shock temperatures, calculated con- 
tinuum temperatures release paths, and the calcite devolatilization equilibria 
are shown in Figure 14. This figure shows the P-T phase diagram for CaC0 3 , the 
Hugoniot (continuum) P-T data for the Solenhofen limestone and Dover chalk 
reported here, and the measured shock temperatures of single crystal calcite of 
Schmitt (personal communication) and Kondo and Ahrens [1983]. The calcite VI 
liquid boundary is drawn through the two measured shock temperature values, 
but the conclusions drawn below do not depend critically on this particular 
placement of the boundary. Representative isentropic release paths, using the 
parameters in Table 3, are shown for several points. The diagram illustrates 
that the isentropic release paths of shocked non-porous carbonates (Solenhofen 
limestone), calculated using continuum temperatures, do not cross into the 
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CaO(,) + C0 8(Ttpor ) field until very low pressures (10~ 7 - 10 _fl GPa) are reached. 
Shocked Dover chalk, particularly at shock pressures as high as 10 GPa, relaxes 
into the solid plus vapor region at pressures of 10 _s - 10“ 4 GPa. However, isen- 
tropic release from points corresponding to measured shock temperatures pass 
into the liquid field region at pressures of 10 1 GPa or higher, into the liquid plus 
vapor field at pressures of 10 -1 GPa or higher, and then into the solid plus vapor 
region at pressures of about 10“® GPa. The high measured shock temperatures 
are related to shear banding phenomena, and it thus appears that localization of 
thermal energy is an integral part of the devolatilization process. The diagram 
also indicates that the effect of C0 2 pressure on the fraction devolatilized will be 
more modest in the case of shear band controlled devolatilization than in the 
case of devolatilization controlled by continuum temperatures. 

Conclusions 

Solenhofen limestone shocked into the calcitc VI region retains its shear 
strength in the shocked state at pressures between about 12 and 20 GPa. Equili- 
brium thermodynamic calculation of the energy and entropy required for inci- 
pient devolatilization combined with measured values of the post-shock energy 
and entropy gain fail to quantitatively predict the amount of volatile loss upon 
release. This result indicates that inhomogeneous deformation processes occur 
in this pressure range. Therefore the density of the shear bands created in this 
shock-pressure interval must be insufficient to modify the bulk mechanical pro- 
perties of the polystalline aggregate. The shear bands may be sites of initiation 
of partial melting and nucleation of evolved C0 2 upon release. At shock pres- 
sures above about 20 GPa loss of strength of the material may occur, implying 
that the shear band density is large enough to affect mechanical properties. 

Isentropic release paths of material initially at the shock pressure and the 
measured shock temperature pass through the vapor plus liquid field into the 
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vapor plus solid field, giving further support to the conception that devolatiliza- 
tion is associated with shear band formation. 

At shock pressures of ~10 GPa, Dover chalk (p 0 = 1.40 g/ cm 3 ) appears to 
have lost its shear strength, but at 5 GPa the results are ambiguous. Post-shock 
entropy calculations indicate that greater than 905 of the C0 2 is devolatilized 
upon release from 10 GPa pressure. 

Equilibrium calculations and experimental observations indicate that the 
amount of devolatilization of volatile-bearing minerals upon impact is dependent 
on the ambient partial pressure of the volatile species. The shock pressure 
required for incipient devolatilization increases with ambient volatile species 
partial pressure. This effect causes the range of planet sizes in which partial 
devolatilization of carbonates occurs to be larger than that calculated based on 
devolatilization experiments carried out in air [Lange and Ahrens. 1963, 1964]. 
Volatile release is enhanced relative to that calculated by Lange and Ahrens in 
the early stages of planetary accretion, when the atmosphere is sparse, and is 
inhibited in the later stages, when the atmosphere is relatively dense. 
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Table 1. Summary of Intermediate, State, Hugoniot State, 
and Release Path Release Transition Data 



Shock States* 

Release Path Phase 
Transitions 

Shot# Projectile Projectile 
Velocity 

s 

km/ sec 

Initial 

Density 

Pe 

g/cm a 

Shock Particle 
Velocity Velocity 

U « up 

km/sec km/sec 

Pressure 

P 

GPa 

Density 

P 

g/cm* 

Panicle 

Velocity 

“st 

km/sec 

Pressure 

Pi 

GPa 

Density 

P 

g / cm 3 





SakaMa Urn 

— 





-500 

AV > 3 

2.200 . 

2.591 

5.819 

1.487 

22.41 

3.460 

2.44 

6.46 

2.90 



*.017 

£.010 

*.150 

£.018 

£.44 

£.042 

£.10 

±1.27 

£.08 

506 

Fused 

2.158 

2.584 

(4.60) 

(0.030) 

(0.36) 

(2.601) 

1.26 

6.20 

3.11 


quartz 

*.02 

*.009 

£.08 

£.010 

£.12 

£.021 

£.03 

£.60 

Hi 





(4.413) 

(0.40) 

(4.58) 

(2.841) 








*.435 

£.10 

*1.05 

±.079 








5.081 

0.997 

13.10 

3.215 








*.080 

*.015 

£.16 

±021 




ana 

wP 

a * o 3 

2.310 

2.584 

8.073 

1.538 

24.14 

3.481 

2.65 

3.21 

2.87 



*.160 

*.010 

±.120 

£.110 

±1.67 

£.065 

±06 

£.60 

±.04 

600 

AleOa 

1.965 

2.608 

5.741 

1.319 

19.75 

3.366 

2.11 

4.97 

2.95 



*.005 

*.002 

£.008 

£.004 

±.06 

£.004 

£.04 

£.51 

£.03 

602 

Fused 

1.992 

2.604 

(5.706) 

(0.030) 

(0.45) 

(2.618) 

1.34 

5.97 

3.00 


quartz 

*.010 

.005 

*.300 

£.010 

£.15 

£.021 

£.04 

£.47 

£.02 





(4.774) 

(0.200) 

(2.56) 

(2.715) 








£.09 

±.02 

±30 

£.025 








4.903 

0.936 

11.96 

3.219 








£.000 

£.008 

£.09 

£.013 








Dover Chalk 





! 509 

A 1,0 3 

2.078 

1.434 

3.750 

1.745 

9.39 

2.683 






*015 

*.030 

£.100 

£.018 

£.26 

±.085 




501 

A 1 8 0 3 

2.331 

1 434 

3.845 

1.959 

10.80 

2.924 






±.020 

*.030 

*050 

±.020 

£.22 

£.066 




620 

Alf 0 a 

1.619 

1.365 

2.833 

1.414 

5.66 

2.835 






*.019 

£.024 

£.047 

*0.16 

£.13 

£.064 




622 

ig 0 3 

2.235 

1.365 

3.652 

1.900 

9.47 

2.646 



1 



£.022 

£.024 

£.127 


£.32 

£.131 





♦Higher, pressure rale is final shock rate. All others are intermediate rates 


Values in parentheses are less reliable data. 











Table 2. Summary of Polycrystalline Calcite Hugoniot Data 


Material 

■Reference 

p*. g/cm* 

C a , km/sec 

s 

Pressure Range, 
GP* 

Solenhofen limestone 

this work 

2.504 

3.269 

1.796 

12-24 

Dover chalk 

this work 

1.40 

.67 

1.60 

5-11 

Polycryst ftilme calcite 

Adadurovetal., 1061 

2.703 

3.40 

2.00 

5-13 




3.90 

1.32 

13-51 

Polycrystaliine calcite 

Kalashnikov et al., 1973 

2865 

3.70 

1.44 

10-94 



2.020 

1.74 

1.61 

13-71 



1.706 

1.15 

1.60 

10-59 

Solenhofen limestone 

van Thiel et al., 1977 

2.565 

3.62 

1.39 

6-90 
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Table 3. Properties of the Calcite Polymorphs 



I 

n 

in 

VI 

Zero-pressure density (300 K) p* g/ cm* 

2.71* 

2.71* 

2.71* 

3.0-3. l b 

Bulk modulus K, GPa 

71.1* 

14.7* 

51.7* 

75-95 b 

dK/dP * K' 

4. 15* 

4.62* 

8.28* 

4.1-3.5 b 

Transition energy Etr, J /g 

- 

- 

- 

200- 20 b 

Gruneisen parameter 7 

1.5 (V/V 0 )* 

1:5 (V/V 0 ) e 

1.5 (V/V 0 )* 

1.5 (V/V 0 )* 

Poisson's ratio a 

0.25° 

0.25* 

0.25* 

0.25* 

K„ d in Vp/dP = D p 

2.13* 

0.44* 

1.55* 

1.0-1. 6 d 

Ko. d in V,/ dP = D, 

2. 13* 

0.44* 

1.55* 

1.0-0. l d 

dVp/dT, km /sec- K 

-3xl0 _4d 

-3X10 -44 

-3xl0 -4d 

-S-OxlO -44 

dVg/dT, km/sec-K 

-3xl0 -4d 

-3xl0 _4d 

-3X10" 4 d 

-2.4X10 -4 d 

Specific heat capacity C,, J/g-K 

1.25* 

1.25* 

1.25* 

1.25* 


* Singh and Kennedy, 1974 
b Vizgirda and Ahrens, 1982 

* assumed 

d Anderson et aU . , 1968; Sammisetol., 1970 

* Peselnick and Wilson, 1968; Wang and Meltzer, 1973 

'■ Duiong-Petit value, T * 650; Vizgirda and Ahrens, 1982 



Table 4. Calculated continuum Hugoniot temperatures Tr. 
extrapolated complete release density p e0 , and post-shock 
energy increase AEp 


Shot 

T«.K 

P 00 - g/ cm 8 

AE p , erg/ gxlO 



Solenhofen limestone 


602 

551 * 82 

2.6B±.02 

1.11* .10 

596 

848± 31 

2.81±.03 

2.51* .16 

600 

825 ± 94 

2.68*. 03 

2.40± .06 

590 

945*100 

2.59*. 10 

2.15± .63 

598 

1024 * 98 

2.?7±.03 

4.46* .31 



Dover Chalk 


620 

1089* 40 

2.44±.02 

9.13* .50 

622 

1747± 98 

2.55±.03 

16.21*1.26 

801 

1832* 97 

— 

— 



Figure Captions 


Figure 1. Schematic diagram of particle velocity gauge experiment. Top view, 
active element of foil gauge is oriented in and out of plane of page. Gauges are 
epoxied between plates of Solenhofen limestone and their electrical leads each 
go to an oscilloscope. Impact-induced velocity of gauge generates a voltage 
which is proportional to the magnetic field (generated by the Helmholtz coils). 
Pin signal triggers oscilloscopes. 

Figure 2 . Chalk (porous calcite) target assembly. Electrical leads from gauges 
and pins are omitted for clarity. 

Figure 3. Schematic diagram of current-generating and signal-measuring cir- 
cuits. Only one of four signal circuits are shown. 

Figure 4. Particle velocity versus time after impact f .n shot 600 obtained from 
digitized voltage versus time oscilloscope records. Gauge 1 is nearest to 
impacted surface, gauge 4 is at free surface. Shock velocity was 5.741 km/sec, 
Hugoniot particle velocity was 1.319 km/sec, and shock pressure was 19.75 GPa. 

Figure 5. Position-time (x-t) diagram of a particle velocity experiment. Projec- 
tile approaches from left and impacts stationary target at x=0, t=0. The 
diagram shows the (undesirable) case in which gauge 1 is intercepted by the 
rarefaction propagating forward from the projectile rear surface before arrival 
of rarefaction from rear surface. 

Figure 6. Shock wave velocity U, versus particle velocity u p data for experi- 
ments reported here and comparison with previous results. Filled symbols show 


present data for the Solenhofen limestone, open symbols show data for Dover 
chalk. Crosses with no symbol represent intermediate states. Lines represent 
fits of form U, s Co + s Up. Parameters for Cits are listed in Table 2. Cur ves are 
labelled as follows: 1-Solenhofen limestone, this study, p 0 = 2.594; 2-single cry- 
stal aragonite, [Vizgirda and Ahrens, 1982], p 0 = 2.93; 3-polycrystalline calcite, 
[Adadurov et ai., 1961], p 0 = 2.703; 4-Solenhofen limestone, [van Thiel et a l., 
1977], p 0 = 2.585; 5-polycrystalline calcite [Kalashnikov et ai., 1973], p 0 = 2.885; 
6-polycrystalline calcite, [Kalashnikov et at., 1973], p 0 = 2.02; 7-polycrystalline 
calcite, [Kalashnikov et ai., 1973], p 0 = 1.705; 8-Dover chalk, this study, p 0 = 
1.40. Units of density, g/cm 3 . 

Figure 7. Hugoniot states and release paths for Solenhofen limestone and Dover 
chalk. Filled symbols are for Solenhofen limestone, open symbols are for Dover 
chalk, boxes represent initial densities. Heavy solid lines are fits to data (Table 
2). 

Figure 8. Particle velocity versus time after impact for shot 596, showing multi- 
ple shock fronts on leading edge of record (including HEL), and indicating pres- 
ence of an intermediate release state during unloading. 

Figure 9. Eulerian sound speeds, Cg, on the Hugoniot of shocked Solenhofen 
limestone. Ranges of compressional wave velocity V p and bulk sound speed V p 
for calcite VI are computed using Eulerian finite strain theory (see text) and 
parameters listed in Table 3. Values plotted are for an isentrope originating at 
750 K at zero pressure. Dashed line is along the Hugoniot of calcite VI calcu- 
lated according to equation (18). 



Figure 10. Eulerian sound speeds C E on Hugoniot of Dover chalk. Solid and 
dashed lines represent compressional wave velocities and bulk sound velocities, 
respectively, of calcite polymorphs along a 750 K isentrope calculated using the 
parameters in Table 3. Roman numerals refer to the assumed calcite 
polymorph. Dotted line labelled H represents bulk sound speed along low pres- 
sure (5-13 GPa) polycrystalline calcite Hugoniot of Adadurov at a l., [1981] (Table 
2) using equation (18). 

Figure 11. Post-shock energy increase AE p as a function of shock pressure. Also 
shown is energy of incipient vaporization Ejy calculated for several values of P COg . 

’STP air’ refers to PcOg in dry air, 3.3xl0 -4 bars; 'Tank' refers to P C o g in impact 
tank under shot conditions, about 4.3 x 10“ e bars. Hatched regions represent 
values of AE p calculated using calcite VI parameters (Table 3) for Solenhofen 
limestone and calcite I and II parameters for Dover chalk. The Hugoniot pres- 
sure of the porous material was calculated using equation (27). 

Figure 12. Shock-induced C0 2 loss as a function of shock pressure in calcite as 
determined in experimentally shock loaded single crystal calcite in an ambient 
atmosphere of air. Data are from Lange and Ahrens [1984]. Also shown are C0 2 
losses calculated using shock entropy method for various values of ambient C0 2 
pressure. 'Tank' refers to the ambient C0 2 pressure in impact tank, about 4.3 x 
10 _a bars; 'STP air' refers to the P^ in dry air at standard temperature and 

pressure conditions, about 3.3 x 10 -4 bars; Solid lines refer to the following com- 
bination of parameters of calcite VI; p 0 = 3.1 g/ cm 3 , = 95 GPa, K' = 3.5, 
E lrtni = 20 J/g. For dashed lines p 0 = 3.0 g/ cm 3 , K p = 75 GPa, K’ = 4.1, E trani = 
200 J/g 


Figure 13. Shock entropy as a function of shock pressure. The curves for 
porous calcite are calculated using the calcite Hugoniot of Adadurov at a l. 
[1981], the correction for porosity embodied in equation (27), and the equation 
of state parameters for calcite VI listed in Table 3. Scv and Sjy are calculated as 
a function of Pco, as described in text. 

Figure 14. CaC0 3 phase diagram showing schematic paths of isentropic 
decompression from shocked states. The P*T phase diagram is based on the 
results of Carlson [1983], Huang and Wyllie [1978], Irving and Wyllie [1973], 
Baker [1982], and Jamieson [1957]. Stability fields of each polymorph are 
labelled with the appropriate roman numeral. The range of Hugoniot and 
release path temperatures for Solenhofen limestone are calculated using the 
range of values of properties of calcite VI (Table 3). Hugoniot and representative 
release path temperatures of Dover chalk are based on the low pressure calcite 
properties (Table 3). Open squares are shock temperature measurements on 
single crystal calcite by Kondo and Ahrens [1983] (K & A) and by Schmitt (per- 
sonal communication. 1984; S). The representative isentropic release path for 
the point S is calculated using the calcite VI physical property data. Note that 
10 -4 GPa corresponds to 1 bar pressure. is the invariant point, at which 
CaC0 3 ( g ), CaO(,), C0 2 ( vtpol .), and liquid coexist. CaC0 3 melts incongruently between 
points Q, and Q 2 . and melts congruently at pressures higher than Qg. Metastable 
calcite !I and III fields are separated by dashed lines. Extrapolated or inferred 
phase boundaries are shown by dash-dot lines. Dotted line represents the ara- 
gonite 1 calcite VI phase boundary constrained to lie below the limestone 
Hugoniot. Error bar on the dotted line represents the position of the line calcu- 
lated using the range of property values for calcite VI listed in Table 3. 
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